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NEUTRAL AIR WAVES IN THE THERMOSPHERE 

by 

H. Volland 
ABSTRACT 

Under the influence of gravity, heat conduction, molecular viscosity, Coriolis 
force and ion drag eight plane characteristic waves obliquely incident on a 
horizontally stratified atmosphere can propagate, four of them upward and the 
other four downward. The four pairs of characteristic waves are the well known 
acoustic-gravity waves, heat conduction waves, and the ordinary and extraordinary 
viscosity waves. The eigenvalues of the characteristic waves are calculated nu- 
merically giving the following main results: 

(a) Under the influence of heat conduction acoustic-gravity waves exist in 
the whole frequency range and at all angles of incidence. 

(b) Coriolis force and ion drag make the atmosphere anisotropic with respect 
to the characteristic waves. Their propagation characteristics at east- 
to-west and at north-to-south propagation differ from each other. 

(c) At the equator upgoing internal gravity waves prefer a predominant 
west to east path while downgoing heat conduction waves tend to propa- 
gate from north to south. 

(d) Downgoing internal gravity waves travelling from north to south lose 
their abnormal behavior under the influence of the Coriolis force and 
have phase and group velocities the vertical component of which go 

in the same direction. 
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NEUTRAL AIR WAVES IN THE THERMOSPHERE 


I. INTRODUCTION 

It is a well known experience that the air within the lower atmosphere be- 
haves as a heat insulator and is free of molecular friction. The negligibly small 
heat flow within air is responsible for the adiabatic behavior of acoustic waves 
and leads to an "adiabatic" velocity of sound 


C a = vyfi T 


( 1 ) 


in which y = c p /c v is the ratio between the specific heats at constant pressure 
and at constant volume, /x =R /M is the ratio between the ideal gas constant R 
and the molecular weight M, and T is the mean temperature of air. 

In order to estimate the velocity of heat flow within quiet air we use the one 
dimensional equation of heat conduction 


p c 

“ V 



3 2 T 
3 z 2 


( 2 ) 


in which ~j5 is the mean density, k is the coefficient of heat conductivity, t is the 
time and z is the altitude coordinate. Replacing t by a characteristic time r 
and z by a characteristic length for which we choose the scale height H, equation 
(2) leads to a velocity of heat transport 


V 


a 


H_ = K 
t ~Pc v H ‘ 


( 3 ) 


Using the numerical values at the earth's surface 



* 


K = 1.8 X 10 2 = 3.0 x 10 3 erg/cm sec °K 

p = 1.29 x 10“ 3 g/cm 3 
c v = 7.2 x 10 6 erg/g °K 


H = 8 km 


gives a velocity of 


V g = 4 x 10“ 7 cm/sec 


verifying the adiabatic behavior of air at low heights. Since the coefficient of 
heat conductivity k only depends on the temperature T the situation reverses 
completely if we consider the upper atmosphere. At a height of 200 km mean 
values are 


T = 1000 °K 
P - 2.7 x 10“ 13 g/cm 3 


H = 40 km 


yielding from equation (3) a velocity of 


V g = 8 m/s ec 


which is by no means negligibly small. 

The coefficient of molecular viscosity v is proportional to the coefficient 
of heat conductivity k and therefore also depends only on the temperature. The 
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viscosity usually is characterized by the Reynolds number which is defined as 


= l L 2 
r 77 


where L and r are again characteristic length and time. Replacing L by a 
reciprocal wave number C a /co and r by the period l/a> of a wave with angular 
frequency c*> and a phase velocity equal to C we obtain from equation (1) 


r ,Z£. 

a COT) 


( 4 ) 


Here we used the ideal gas equation 


P = fJ-P T 


( 5 ) 


in order to introduce the mean pressure p . 

By choosing the same numerical values as above, taking (Chapman and 
Cowling, 1959 ) 


V = 


2.5 c v 


and considering the angular frequency o> = 1 sec * 1 we obtain 


8.5 x 10 9 at sea level 


R = 

a 


3.8 at 200 km height . 
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Thus with respect to the influence of viscosity we again find an entirely different 
behavior between the lower and the upper atmosphere. The range of turbulence 
is related to high Reynolds numbers and small Richardson numbers. In heights 
below 100 km wind shears and thermal instabilities due to the small heat con- 
duction are responsible for turbulent mixing of the different constituents of the 
atmosphere. In heights above 120 km the air is nearly in diffusive equilibrium 
as the result of negligibly small turbulence. 

In the following paper we consider plane atmospheric waves propagating 
obliquely through a horizontally stratified atmosphere under the influence of 
heat conduction, molecular viscosity, ion drag and the Coriolis force. We shall 
see that within the upper atmosphere three important parameters exist closely 
related to the parameters C a , V a and R a of equations (1), (3) and (4) which are 


C = 4t “isothermal" velocity of sound 


V = heat conduction velocity 

, c / H 


( 6 ) 


R = _H_ Reynolds number 

CO 7} 


Within the lower atmosphere the two well known acoustic-gravity wave modes 
(one upgoing and one downgoing) exist (Hines, 1960). Their typical phase velocity 
is C in the high frequency region. Under the influence of heat conduction two 
new modes are created, called heat conduction waves (Volland, 1967). Their 
typical phase velocity is V in the low frequency region. Molecular viscosity 
generates four additional modes, ordinary and extraordinary viscosity waves. 
Their phase velocity is typically C/^R. 

In section 2 the eigenvalue equation is derived. Section 3 deals with special 
analytical solutions within an isothermal atmosphere. The difference in the be- 
havior of the wave modes as compared with the original adiabatic acoustic- 
gravity waves introduced by Hines (1960) is stressed. In section 4 the general 
dispersion equation of all 8 modes is discussed in greater detail. 

The first numerical calculation of vertical wave motions in the upper at- 
mosphere taking into account heat conduction was made by Harris and 
Priester (1962) in their famous Harris-Priester-model. Pitteway and Hines 
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(1963) and Midgley and Liemohn (1966) discussed the influence of heat conduction 
and molecular viscosity on the acoustic-gravity waves. In this paper the prop- 
erties of all kinds of wave modes within the upper atmosphere are examined. 


II. THE EIGENVALUE EQUATION 

We start from the equations of conservation of mass, momentum and energy 
and the ideal gas equation of the neutral gas which are 


Bt 


+ div (p v) = 0 


dv 


p — + div o- + v p (v - v.) -20.x v + grad p - p g = 0 
dt 1 


(7) 


J np 

p -j-i- + p div v + /3 - div ( k grad T) = 0 


p - p. p T = 0. 


p density; v = (u,v,w) velocity 
P pressure; T temperature 

v number of collisions between neutral molecules and ions 
o ■ viscous stress tensor 

f Bvj Bv k 2 -.1 

a.. = - 77 < + o\ , div v r el ement of cr 

lk '{ Bx k Bx. 3 lk J 

77 coefficient of molecular viscosity 
0 earth's rotational vector 
g gravitational acceleration force 
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e v specific heat at constant volume 


- <v • V B 0 


ion velocity, supposed to be zero when orthogonal to the geomagnetic field B 0 



v. 


Bx 


k 


viscosity heating 

k coefficient of heat conductivity 
/x = R/M 
R gas constant 
M molecular weight 

For solving the system of equations (7) a consequent perturbation method is 
applied assuming that the time independent mean values like density p, pres- 
sure p and temperature T are already known and that the mean velocity of the 
air v is zero. The perturbation is considered to be a plane harmonic wave of 
angular frequency co and of wave number k obliquely incident on a horizontally 
stratified plane atmosphere in which all parameters depend only on the vertical 
component z . Then all variables are functions of the coordinates x , y , z , and 
time t according to 


f (z) e j&,t ~ i k sin 6 


(x cos A + y sin A). 


(8) 


6 is the angle of incidence and A is the azimuth of the plane of incidence of the 
wave with respect to geographic south. We assume that no internal energy 
source exists. Therefore, the vertical component of energy transported by the 
different wave modes is continuous at any internal boundary. This implies 
Snell's law: 


6 



I 


k sin 6 - k Q sin = const. 
A = A Q = const. 


( 9 ) 


j 

It is possible, therefore, to relate wavenumber k, angle of incidence <9 and azi- 
I muth A of a wave to an arbitrary wave with wavenumber k Q =a >/C 0 and angle of 

incidence d Q at the arbitrary height z 0 travelling in the A,, direction. We choose 
as such an arbitrary wave an undissipated acoustic wave which has in the frequency 
limit co — ix within a loss free atmosphere the isothermal acoustic phase velocity 
C 0 = / m T 0 of equation (6) and the real wave number k 0 . The angle of incidence of 
this wave is real. C 0 plays the role of a reference velocity equivalent to the ve- 
locity of light in the electromagnetic wave propagation in plasma. 

For convenience we rotate the coordinate system in such a manner that the 
new x-axis points in the A 0 direction. Then the horizontal components of the 
velocity with respect to the new coordinate system are 


u = u cos A Q + v sin A Q 


'V 

v = -u sxn A 0 + v cos A 0 . 

i 


The values u, v, w, Ap, A p and AT are the deviations from the mean values and 
are considered to be small compared withp (z), p(z) and T (z) of the quiet at- 
mosphere. Thus all products of these small values can be neglected. We nor- 
malize these variables according to 



( 10 ) 


The derivation with respect to z has been replaced by a prime (9/Bz = '). 


! 
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Introducing the expressions of equation (10) into the system of equations (7) 
and considering the time and spatial dependence of equation (8) the variables w" 
and A p can be eliminated which leads to a system of first order differential 
equations which in concise matrix form is 


e' - jk Q Ke = 0, (11)’ 


where 


e (z) = 



is the row matrix containing the 8 independent variables defined in equation (10). 
The coefficient matrix 



has the submatrices 


8 



0 


Here 


-J 


K 4 = 


Rdj (1 - j B ? ) - j Sq 0 


jRB 6 d 1 0 


-) 


-jRB s d 1 0 Rdj (1 - jB g ) - jSg 0 


the following abbreviations have been used: 


k - — 

0 C 
'“0 


C 0 =4t>,„ 


S 0 = sin e o 


C = vV T(z) 


A = 


2k 0 H 


H =^-L= - — 


1 2 k oHl 

1 


A ’ = 2k 0 H 2 


h --F 


T 

H 2 = r=r~ 

T' 


G = 


2AC„ 


V = 


SPH 


d ! = 


c 1 2 

F 


j _ 7 

2 " (7-1) 


c 

v 


R =- 


CO T) 
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1 


s = 




y _ 2n 7 _V 
1 ~ 2 - — 


•a i / 3 2A! 

D = 1 + |RZ 2 cos 2 I + i ^ f S 2 - 4A 2 - -i 


B, 


B„ 


= Zj sin iJ sin A Q ± Z 2 sin I cos I cos A 0 


B, 


B, 


= Zj s in & cos Aq ± Z 2 s in I cos I s in A Q 


B, 


B, 


= Zj cos d ± Z 2 cos 2 I cos Ap sin Ag 


B, 


B„ 


cos 2 Ajj 


= Z 2 1 1 - cos 2 I 


SI 


in 2 Aq 


T sin iJ . T 2 cost? 

cos I = sin I = 


/l + 3 cos 2 t? /l +3 cos 2 t? 


t? geographical co-latitude 
I geomagnetic dipangle 
A 0 azimuth 

In this calculation the geomagnetic field has been approximated by a dipole 
field with its axis parallel to the earth's rotational axis. Moreover, the deriva- 
tives of fj . , k and 7] with respect to z have been neglected. 

In order to solve equation (11) uniquely it is necessary to know the conditions 
at the upper or lower boundary of the model atmosphere. Physically appropri- 
ate solutions need the separation between upgoing and downgoing waves at these 
boundaries. Such waves are the characteristic waves which are defined by the 
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eigen values q v of the matrix K . They can be found from the eigenvalue equa- 
tion 


|K-\„ E| = 0, 


(13) 


with 


= -(q + jA) 


and 


E unit matrix. 


Now a characteristic wave has the z-dependence 


f „( z ) = e 


j k 0 X v 2 


= e 


• jk 0% z + 2H 


(14) 


Here the real exponential factor 1/2H = k 0 A has been split from the eigenvalue 
q y merely because of convenience. It gives the increase of the amplitude with 
height of a non-dissipative wave in an isothermal atmosphere. Therefore the 
imaginary term of q v is negative or zero for an upgoing wave. 

A characteristic wave is uncoupled from the other waves if the elements of 
the matrix K and therefore the q v are height independent. Since K depends on 
V and R which themselves depend on ~p and p , respectively, the characteristic 
waves are always coupled in a realistic atmosphere. Therefore, in order to find 
out the relative importance of the different wave modes one has to solve the sys- 
tem of equations (7) by numerical methods beginning at an arbitrary level z 0 and 
integrating upwards (or downwards) to a second arbitrary level z n (see figure 1). 
This leads to a solution of the general form 


e ( Z n> = T Z " e < Z 0> 


(15) 
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The physical quantities e (z 0 ) and e (z n ) at the boundaries are connected with 
the characteristic waves by the equation 

e = P c . (16) 


Here 



is composed of four upgoing waves a u and four downgoing waves b v while P 
has to be chosen in such a manner that 



P* 1 K P = N = 




V 


\° • • • ° \! 


(17) 


is the diagonalized normal form of the matrix K. 

The elements of the eigen vectors P„ (P is the v th row matrix of P ) are 
determined from the system of linear equations 


(K-\„ E) • P v = 0. 


(18) 
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Finally by well known matrix calculus (Volland, 1968) one obtains from 
equations (15) and (16) the scattering matrix M which is defined by the rela- 
tion (see figure 1) 



The scattering matrix completely describes the behavior of the atmosphere 
with respect to plane atmospheric waves. 

If the coupling between the different modes is weak a "ray optics" approxi- 
mation is valid in which each mode is considered to propagate without coupling 
with a height dependence according to 


" jk o f z z 0 * (?) d * + /,* 2i^> 


(19) 


In the complex eigenvalue 


% = %-}P v (20) 

the imaginary part /3 is responsible for the attenuation of the wave mode in the 
vertical direction while the real part a p is a measure of the phase velocity V p 
of the wave: 






n 


( 21 ) 


where n is a real refractive index. 

The direction of the wave normal is given by 


d p | = arctg 

V 



( 22 ) 
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while the ray direction can be found from the relation 


tg d. 


Xg 


3a sin A. 3 a 

l v 0 v 

cos A- + 

0 3S„ c 


3 A„ 


tg^ygl =- sinA 0 


3 a cos A„ 3 a 

v o v 


3S„ 


3 A„ 


( 23 ) 


where i? xg and d yg are the angles of inclination of the ray components in the 
(x,z)— and the (y,z)-plane with respect to the z-axis. 

The ray propagates with vertical group velocity 


V | = E° (24) 

g v 3(aa v ) 3a v 

3a> 0 3 S q 

thus allowing "ray treatment" of the waves in complete analogy to electro- 
magnetic wave propagation (see e.g. Budden, 1961). The WKB-approximation of 
ray optics is distinguished by a specific normalization of the matrix P in such a 
manner that 


/ 


0 • • 


(P- 1 P') = 



• • 0 


(25) 


is a matrix in which the diagonal line consists only of zeros. (See e.g. Volland, 
1968). 

The time averaged vertical component of the energy transported by the v th 
wave mode is approximately (see e.g., Eckart, 1960) 


S z I =i-Re (pw) | = ip 0 w 0 (g 

V ^ V ** 


lv ^2^ ® \v 


a a* 

V V 


b b* 

V V 


(26) 
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where e. and a are the elements of the two first lines of P and the star in- 
dicates complex conjugate values. In the case of the WKB-approximation the 
amplitudes a v a* and b v b* of the different modes are directly comparable. 

m. SPECIAL SOLUTIONS OF THE EIGENVALUE EQUATION 
(ISOTHERMAL ATMOSPHERE) 

Numerical calculations showing the relative importance of the different wave 
modes shall be given in an additional paper. Here we confine ourselves to the 
discussion of the eigenvalue (or dispersion) equation. In this section a few 
special cases are considered which lead to analytic solutions of the eigenvalue 
equation. Such solutions are essential because they allow one to identify the vari- 
ous wave modes in more general cases treated numerically in the subsequent sec- 
tion. For convenience we consider only an isothermal atmosphere. Then 
T' = 0 and therefore A 2 =0, A x = A and = 1 in K (equation (11)). 

In the case of a wave propagating at the equator (d = 90°) in an east- west 
direction (A 0 = ±90°) 


B 3 - B 4 “ B 5 “ B 6 - B 8 - 0 


and 


B, 


B = ± Z. 


B, 


Then the eigenvalue equation (13) has the form 


= 0 


( 27 ) 
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CO | K> 


with 




(29) 


Now the eigenvalue equation can be separated into 

|K S | |KJ =0 


(30) 
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which implies 


|K S I =0 

or (31) 

IK, I =<> 


The equation | K g | =0 gives the eigenvalues of a transversal type of viscosity 
wave which we name the extraordinary viscosity wave mode: 


» v 2 ti R (i-i|)=0 


q 7 = T j + jR -jA 
8 

The upper sign stands for the upgoing wave, the lower sign for the downgoing wave. 
This wave is a pure transversal wave inasmuch as only the horizontal component 
v orthogonal to the direction of wave propagation is involved. The wave is heavily 
attenuated for large values of R and is a purely evanescent one at R = 0. It is 
completely uncoupled from the other wave modes in this direction of propagation. 
At vertical incidence its phase velocity is (equation (21)). 


/2C 0 

y =± - 

P 


- + 



and its group velocity is (equation (24)). 


The remaining equation | K s | =0 can be solved analytically for two special 
cases. 
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3(a) Small Reynolds numbers (R - 0, S - 1) 


The equation |K S | = 0 can be written as 


(\ 2 + S 2 ) 2 jx(\ + 2jA) + S 2 + i£j= 0 (33) 

and has the solutions 


Qi = Q s = + i s o - i A 
2 6 

(34) 

q 3 = 4 j ]/a 2 + S 2 + ± T j A (l +— — ) 

4 0 7 2yA \ 2 A 2 / 

where the last approximation is valid in the low frequency region (A 2 » G » 1). 
The upper sign stands for the upward going waves, the lower sign for the down- 
ward going waves. The acoustic -gravity mode (index 1 and 2) and the ordinary 
viscosity mode (index 5 and 6) are purely evanescent waves while the heat con- 
duction mode (index 3 and 4) has the phase velocity (see equation (21)) 


V = x — 

]/ t vis .y 

in the low frequency region which in the case of vertical incidence (S 0 =0) is 
equal to the phase and group velocity (see equation (3)) 


v a =yv. 
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Thus the heat conduction wave mode is the only carrier of energy. 


The ray direction in the low frequency region depends only slightly on S Q . 
According to equation (23) the heat conduction mode therefore tends to propagate 
purely vertically. The amplification factor of the upward going wave is negative 
and small compared with the factor A in equation (14): 


-fi 3 + A ~- 



showing that part of the energy of the upward going heat conduction wave is con- 
verted into internal energy of the surrounding air. The same is true for the 
downward going component which is attenuated according to 


-/ 3 4 + A ~ 2A 


a factor twice as large as the attenuation factor of a non dissipative downgoing 
wave (see equation (16)). The effective attenuation factors of the other wave 
modes according to equation (32) and (34) do not depend on the earth's gravita- 
tional field: 


Im = Pv - A - ± S 0 


O = 1,2, 5, 6, 7,8) 


3(b) Large Reynolds numbers (R - oo, S 0, RS - 4/3 j) 


In order to solve this problem we have to transform the matrix K g in such a 
manner that we can split the determinant |K S | into two independent determinants: 


K 


s 


Kj + k E 0 

o k 7 
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We do this by multiplying the 5th row of K g (equation (28)) with 


± 



and 


S 


o 


1 


i Z 2 


respectively, and add this row to the 1st and 2nd row. Then the off diagonal sub- 
matrices of K s no longer depend on R p and their elements become negligibly 
small compared with the elements of the matrix K ? when R goes to infinity. 

Thus the eigenvalue equation becomes 


|k 7 | |k x +^e| = o 


and we find the eigenvalues of the ordinary viscosity mode from 



= 0 as 


q 5 = T j i S 2 + RZ 2 + jR- jA (35) 


while the solution of |K X + A.E| = 0 is 


q x = T /a- 7b 

(36) 


q 3 = T + /b 

4 


with 


a=|(F- jG)-A 2 -S 2 


b=l(F- jG) 2 + ^{f - (7 - 1) 


2AMS 0 '[ 

U J 
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M 


F 


= ±Z X - 2AS 0 


= U + 


Z X M 

“ 0 “ 


G 


C 2 

wVH 



U = 1 - j z 2 


s 0 = Sin e 0 


z 1 = 


2Q 

O) 


z 2 = 


O) 


In figure 2 the real parts a x and a 2 of the eigenvalues of the upgoing acoustic- 
gravity wave and of the heat conduction wave calculated from equation (36) have 
been plotted versus the sine of the angle of incidence S 0 = sin 8 0 for different 
angular frequencies between co = 10 -5 and 10 2 sec - and for neglected Coriolis 
force and ion drag (Z t =Z 2 = 0). The Harris-Priester model 5 at 12 00 local 
time and 200 km height (CIRA, 1965) has been chosen as the atmospheric model. 
For convenience the ordinate a is given in two different scales separated by a 
dash-dotted line, respectively. Because of equation (36) the eigenvalues of the 
downgoing waves differ only in sign from the eigenvalues of the upgoing waves. 

In figure 2a the characteristic behavior of the acoustic-gravity waves can 
be seen which in the low frequency range (co < 10 -2 sec -1 ) has a downward 
direction of phase propagation while the energy of the wave propagates upward. 
This behavior can be derived from equations (22) and (23) in which the angles of 
the direction of the wave normal and of the ray (in our special case the ray 
direction does not depend on A 0 ) are given as: 




S 


o 


a 


tgtf g 


~d<x 


3S 


o 


Figure 3 a shows the geometry of the vectors of the wave normal and of the 
ray direction derived from equations (22) and (23) at two frequencies belonging 
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to an acoustic wave type (co = 10“ 2 sec - and to an interval gravity wave type 
(a; = 10“ 4 sec -1 ) and illustrates again the reversal of vertical phase and group di- 
rection of the internal gravity waves. But contrary to the adiabatic case of zero 
heat conductivity (see Hines' (1960) fig. 9) ,the interval gravity waves as well as 
the acoustic waves, exist within the whole range of S 0 . For comparison the 
adiabatic waves ( k = 0) at the same frequencies have been drawn in fig. 3a as 
dashed lines. There exists a continuous transition from the interval gravity 
wave type to the acoustic wave type again in contrast to the adiabatic case where 
no wave mode can exist within the frequency range of 



At high frequencies the acoustic wave approaches the phase velocity of sound 


The heat conduction modes in fig. 2b only slightly depend on the angle of in- 
cidence. These wave modes therefore transport their energy nearly vertically 
(see fig. 3b). In the low frequency range (co < 10 -4 sec ~ ) the phase velocity of 
the heat conduction wave is (see equations (21) and (6)) 


V 


p 





V 


a 


and approaches the heat conduction velocity V . In the high frequency range 
the heat conduction modes become evanescent. 


The dependence of a l and a 3 (and therefore of a 2 and a 4 ) on the parameter 
of the upper atmosphere at vertical incidence (S 0 = 0) has already been shown 
in a previous paper (Volland, 1967). There are some slight differences in the 
denoting of the variables in that paper compared with the present paper: 


old notation new 
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Moreover the following printed errors should be corrected: In fig. la ( z = 400 
km) of (Volland, 1967) the trace of n 4 should be moved upward by a factor of 
10, and the element A 21 of A in equation (39) should be read as 


-jcoa 21 (1 - SD). 


Figure 4 gives the imaginary terms k 0 /S 1 and k Q /3 3 of the acoustic-gravity mode 
and of the heat conduction mode. Since according to equation (14) the factor 
1/2H has been split from the imaginary terms of q the upgoing waves have in- 
creasing amplitudes if k 0 /3 < 1/2H and decreasing amplitudes if k 0 /3> l/2H. The 
last one occurs at co > 10 -2 sec -1 for the acoustic-gravity mode. For compari- 
son we note that the adiabatic acoustic -gravity waves have the value /3 1 = 0 over 
the whole frequency range. 

— 3 

The heat conduction mode has an almost constant amplitude at w < 10 
sec -1 while at frequencies co > 10 -1 sec -1 the wave is attenuated nearly com- 
pletely after a distance of some km. 


Under the influence of the Coriolis force and ion drag the atmosphere behaves 
like an anistropic medium with respect to neutral air waves. This can be seen 
from the A 0 -dependence of the ray direction in equation (23). In our special case 
of equation (36) the propagation in the east- west direction is different from the 
propagation in the west-east direction recognized by the change in sign of Z 1 . 

But this difference is rather small. The ion drag tends to decrease the influence 
of the Coriolis force at very low frequencies as can be seen from equation (36) 
where at S Q =0 the factor F has the value 


fiU^) fz v 

''foJ 


Z 2 » 1 


4Q 2 
2 


V. co 


(37) 


Z. »; 


z 2 = 0 


2 -1 
thus decreasing like co~ at Z 2 = 0 as compared with a decrease like co at 

Z 2 ^0. 


Because of our assumption that the ion velocity is zero orthogonal to the 
geomagnetic field lines (see equation (7)) the horizontal ion drag is fully effec- 
tive for the east-west propagation while for the north-south propagation the 
horizontal ion drag completely disappears. Because of equation (37) we there- 
fore expect a rather large difference in the behavior between these two directions 
of propagation. This will be confirmed in the next section. 
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IV. GENERAL BEHAVIOR OF THE EIGENVALUES 

Apart from some special analytical solutions which we gave in the previous 
section, equation (13) has to be treated by numerical methods. A convenient 
program available in the SHARE Program Catalog is SAD 3099 (EIG4) which 
solves eigenvalues of complex matrices and is coded for the IBM 7090/94. This 
program has been used for the following calculations. 

The atmospheric model is again the Harris-Priester model 5 at 12 00 local 
time and 200 km height (CIRA, 1965) which has the data 

H = 43.3 km ; k = 8.64 x 10 3 erg/cm sec °K 

T = 1100°K ; c y = 7.23x 10 6 erg/g °K 

p = 2.71 x 10" 13 g/cm 3 ; 7=1-5 

Moreover the following additional numerical values have been used: 


v = 1 0 -4 sec -1 ; 0 = 90° 

O = 7.23 x 10" 5 sec -1 ; A 0 =9O°, 180° 

m “4 

7? = 4.78 x 10 g/cm sec 


In figures 5 and 6 the calculated real parts of the eigenvalues a 2 , a 4 , a 6 , 
a g of the downward going four wave modes at the equator (t? =90°) have been 
plotted versus the sine of the angle of incidence S 0 = sin 6 0 for the different 
frequencies a > = 10" 5 sec" 1 to 10 2 sec" 1 . The full lines give the modes propa- 
gating from west to east (A 0 = 90°), the dashed lines are due to the propagation 
from south to north ( A 0 = 180°). Again for convenience the ordinates in figure 5 
have different scales divided by the dash-dotted lines. Apart from the acoustic- 
gravity mode the difference in the behavior between the downgoing and the up- 
going waves is rather small. Therefore 
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within less than 10% accuracy, and the a 2 , a 4 , a g and a g are plotted as negative 
values in order to compare them directly with the values a l and a g in figure 3. 

Likewise the east-west propagation and the west -east propagation on the 
one hand and the north-south and the south-north propagation on the other hand 
show nearly the same propagation characteristics. Thus we confine ourselves 
to downgoing waves propagating from west to east (full lines) and from south 
to north (dashed lines). 

Comparing figure 5 with figure 2 we do not see large differences in the be- 
havior of the acoustic-gravity mode at east-west propagation within the low 
frequency range (a> < 10" sec" ) and in the whole frequency range of the heat 
conduction modes. This shows that the Reynolds number R which is rather large 
there (R > 242 for co < 1 0~ 2 sec* 1 ) does not influence very much these 
waves. On the other hand at frequencies u* > 1 sec -1 (R < 0.242) the acoustic - 
gravity mode differs strongly from the R-free type in figure 2 giving evidence 
that in this frequency range the approximation R - 0, treated in section 3(a), be- 
gins to be valid. The acoustic wave becomes an evanescent one there. 

The striking feature of figure 5a is, however, the behavior of the waves 
propagating from south to north (the dashed lines in figure 5). Now the gravity 
modes at co < 10" sec " have lost their abnormal propagation characteristics 
and behave in the normal sense with phase and group velocity going in the same 
vertical direction at angles of incidence S 0 < 0.9. This is essentially due to the 
Coriolis force which is not reduced in its influence by the ion drag at this azimuth 
of propagation (see equation (37)) as was the case for west to east propagation. 

The phase velocity of the heat conduction mode (figure 5b) increases with 
approach to the north-south propagation at low frequencies. In a ray optics 
treatment a heat conduction wave travelling downward will shift therefore its 
direction towards a predominantly north to south direction. This is illustrated in 
figure 7 where the refractive index surface of the upgoing heat conduction wave 


n 
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has been drawn in a (S Q , A 0 )-plane. The vector of the ray direction is perpen- 
dicular to this surface (see equation (23)) and generally is bent out of the plane 
of incidence. The ray will laterally deviate during its propagation till it reaches 
the stable north-south plane where the lateral deviation out of the plane of in- 
cidence is zero. 
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The real parts of the eigenvalues of the ordinary and the extraordinary 
viscosity waves in figure 6 behave very similar to the heat conduction waves 
apart from the lowest frequency range. 

Figures 8 and 9 show the imaginary parts k 0 /3. (the attenuation factors) of 
the downgoing waves at the equator again plotted as negative values in order to 
compare them with the values in figure 4. Here the azimuth of propagation has 
no striking influence on the attenuation factors. Thus only the east-west propa- 
gation is shown. 

While the attenuation factors of the heat conduction waves in figure 8b only 
slightly differ from the viscosity free mode in figure 4b the acoustic-gravity 
waves (figure 8a) are heavily attenuated in the high frequency range (^ > 1 sec” 1 
contrary to the viscosity free case. The reason is that in this frequency range 
R < 1, and the approximation R - 0 treated in section 3 a can be used which 
leads to equation (34) and shows that the acoustic waves become evanescent 
ones. 

The attenuation factors of the ordinary viscosity waves in figure 9 look 
similar to the attenuation factors of the heat conduction waves. But the extra- 
ordinary viscosity waves have even positive values of k 0 /3 g at a> < 10“ 4 sec” 1 . 
This results from the fact that according to equation (32) the effective attenua- 
tion factor of this wave mode is 


-P 7 + aj + A =± yf 

and does not depend on the earth's gravitation field. If A > VR /2 which occurs 
at co = 10” 5 sec” /S 8 becomes negative. 

Figure 10 shows the real part of the eigenvalue of the upgoing acoustic - 
gravity mode at the equator. We see that again west-east and north-south 
propagation differ remarkably in the lowest frequency range though there is 
not such striking reversal in the propagation characteristics as for the down- 
going waves in figure 5. The ray bending as illustrated in figure 7 leads to a 
predominantly west-east propagation for frequencies between co = 10 ” 4 and 10 " 2 
sec" \ 

The same calculations have been repeated for the values d = 45° and 135° 
equivalent to the geographical latitudes of ±45°. The results have essentially 
the same features as figures 5 to 10 but with, of course, different numerical 
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values. A big difference occurs in the acoustic -gravity mode where the reversal 
to a normal behavior of the downgoing wave now appears at the west-east propa- 
gation path while the north-south propagation acts in the usual way. There exists 
no large difference in the propagation characteristics between 45° and -45° lati- 
tude. The equatorial ray deviation of the heat conduction waves toward a north- 
south propagation as illustrated in figure 7 is diminished in ±45° latitude. Like- 
wise at heights greater than 200 km the situation does not change completely 
apart from the fact that with increasing height the Reynolds number decreases. 
Thus the range where the approximation R - 0 (equation 34) is valid shifts to 
smaller frequencies and the heat conduction mode tends to remain the only car- 
rier of energy. 


V. CONCLUDING REMARKS 

It has been shown that eight obliquely incident plane characteristic waves 
can propagate through a horizontally stratified atmosphere, four of them travelling 
upwards and the other four travelling downward. The four pairs of characteristic 
waves are the well known acoustic-gravity waves, the heat conduction waves and 
ordinary and extraordinary viscosity waves. The heat conduction waves and the 
viscosity waves are named according to their relation to finite heat conductivity 
and finite molecular viscosity. Under the influence of the Coriolis force the 
atmosphere behaves like an anistropic medium with respect to atmospheric 
waves. Ion drag tends to decrease the influence of the Coriolis force in such a 
manner that at the equator the west to east propagation path is less effected by 
the Coriolis force than the north to south propagation path. 

Real and imaginary parts of the eigenvalues of the characteristic wave 
modes have been calculated numerically. With the help of these eigenvalues the 
general behavior of the different wave modes has been discussed in detail. The 
main results are the following: 

(a) Under the influence of heat conduction there exists a continuous transi- 
tion from the low frequency range of the acoustic-gravity waves (the 
so called internal gravity waves) to the real acoustic waves at all angles 
of incidence. The internal gravity waves are attenuated transferring 
part of their energy into internal energy of the surrounding air in such 
a manner that the amplitude of an upgoing wave remains nearly constant. 

A new mode — heat conduction waves — is created which shows small 
dispersion within the low frequency range and which transports its 
energy mainly vertically. Group and phase velocity have the same 
direction. 
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(b) Taking into account molecular viscosity two additional modes - ordinary 
and extraordinary viscosity waves - exist. Their propagation charac- 
teristics are similar to the heat conduction waves. 

Molecular viscosity does not change very much the general behavior of 
the internal gravity waves and of the heat conduction waves but alters 
the acoustic waves within the high frequency range into evanescent ones. 

(c) Coriolis force and ion drag make the atmosphere anisotropic with respect 
to the characteristic wave modes. At the equator within the low fre- 
quency range heat conduction "rays" are laterally deviated out of their 
plane of incidence toward a predominately north to south direction while 
internal gravity waves prefer a predominately east to west direction. 
Travelling ionospheric disturbances show some preferred horizontal 
directions of propagation depending on season and latitude (Heisler, 
1963). This might be due to the influence of the Coriolis force in an 
atmosphere varying with geographic latitude and season. If this is true 
it should be possible to decide from the different behavior of the in- 
ternal gravity waves and of the heat conduction waves in a rotating at- 
mosphere whether upgoing acoustic-gravity waves (Hines, 1960) or 
downgoing heat conduction waves (Volland, 1967) are responsible for 
these travelling disturbances. 

Downgoing internal gravity waves travelling north to south lose their 
abnormal behavior under the influence of the Coriolis force and have 
phase and group, velocities the vertical component of which go in the 
same direction. 

The relative importance of the different wave modes for energy trans- 
port remains an open question. Because of the coupling of the different modes 
within a realistic atmosphere this problem can only be solved by a numerical 
calculation using full wave theory. 

In a ray optics treatment one should expect that the time averaged energy 
transport of each wave mode could be calculated from equation (26). But un- 
fortunately equation (26) is only valid under adiabatic conditions. In a dissipa- 
tive atmosphere where heat conduction and molecular viscosity must be taken 
into account the ray direction derived from the wave packet method (equation 
(23)) differs substantially from the direction of the energy flux vector p v (of 
which equation (26) is the vertical component). This situation is equivalent to 
electromagnetic wave propagation where the Poynting vector shows the 
same discrepancy with respect to the ray direction (Hines, 1951). The reason 
for this difference is the fact that the energy flux vector of the atmospheric 


29 


waves contains the gravitational potential (Volland, 1967) which is uncertain by a 
constant factor. The treatment of this problem will be the subject of a forthcom- 
ing paper. 
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FIGURE CAPTIONS 


Figure 1 
Figure 2 


Figure 3 


Figure 4 


Figure 5 


Figure 6 
Figure 7 


. The upgoing (a) and the downgoing (b) characteristic waves at the 
boundaries of an atmospheric model. 

. Real part a of the eigenvalues of the upgoing acoustic -gravity waves 
(figure 2a) and of the heat conduction waves (figure 2b) versus the sine 
of the angle of incidence 6 0 taking into account heat conduction. The pa- 
rameter is the angular frequency co (in sec" ). The dash-dotted lines 
separate different scales. 

. Illustration of the vector of the real refractive index n (which goes 
in the direction of the wave normal) and the vector of ray direction 
(which is perpendicular to the refractive index surface) for an in- 
ternal gravity wave (co= 10" 4 sec -1 ), for an acoustic wave type ( w = 
10" 2 sec -1 ) (figure 3a) and for a heat conduction wave (figure 3b). 

The dashed lines in figure 3a are the equivalent adiabatic acoustic- 
gravity waves. 

. Imaginary part k 0 /3 (attenuation factor) of the eigenvalues of the up- 
going acoustic -gravity waves (figure 4a) and of the heat conduction 
waves (figure 4b). The parameter is the angular frequency co (in sec -1 ). 
Waves with values k 0 /3 ( £yL/2H have decreasing (increasing) amplitudes 
with height. 

. Negative real part -a of the eigenvalues of the downgoing acoustic- 
gravity waves (figure 5a) and of the downgoing heat conduction waves 
(figure 5b) versus the sine of the angle of incidence <9 0 taking into ac- 
count heat conduction, molecular viscosity, Coriolis force and ion drag. 
The parameter is the angular frequency co (in sec -1 ). The propagation 
conditions are: West to east propagation (full lines) and south to 
north propagation (dashed lines) at the equator. The dash-dotted lines 
separate different scales. 

. Negative real part -a of the eigenvalues of the downgoing ordinary 
and extraordinary waves. For details see text of figure 5. 

. Ray refractive index n 3 = j/a^ + of the upgoing heat conduction 
wave in a (S 0 , A 0 )-plane illustrating the bending of the ray out of the 
plane of incidence. 
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Figure 8. 


Figure 9. 
Figure 10. 


Negative imaginary part -k 0 /3 (attenuation factor) of the downgoing 
acoustic-gravity waves (figure 8a) and of the heat conduction waves 
(figure 8b) versus sinus of angle of incidence <9 0 taking into account 
heat conduction, molecular viscosity, Coriolis force and ion drag. 

The parameter is the angular frequency a> (in sec -1 ). The propagation 
conditions are: West to east propagation at the equator. 

Negative imaginary part -k 0 /3 (attenuation factor) of the downgoing 
ordinary and extraordinary waves. For details see text of figure 8. 

Real part a of the eigenvalue of the upgoing acoustic -gravity waves 
versus the sine of the angle of incidence 6 taking into account heat 
conduction, molecular viscosity, Coriolis force and ion drag. The 
parameter is the angular frequency go (in sec -1 ). The propagation 
conditions are: West to east propagation (full lines) and south to 
north propagation (dashed lines) at the equator. The dash-dotted 
lines separate different scales. 
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acousticMgravity wave mode heat conduction wave mode 
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ANGLE OF INCIDENCE sin 8 0 ANGLE 0F INCIDENCE sin0 o 

Figure 2-Real part a of the eigenvalues of the upgoing acoustic-gravity waves (figure 2a) and the heat conduction 
waves (figure 2b) versus the sine of the angle of incidence d Q taking into account heat conduction. The parameter 
is the angular frequency <u (in sec' 1 ). The dash-dotted lines separate different scales. 
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ANGLE OF INCIDENCE sin0 n ANGLE OF INCIDENCE sin 8 




Figure 7-Ray refractive index n 3 = + Sq of the upgoing heat conduction wave in a 

(S n/ A n )-piane illustrating the bending of the ray out of the plane of incidence. 
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ANGLE OF INCIDENCE sin d Q ANGLE OF INCIDENCE sin#. 

Figure 9-Negative imaginary part -k 0 /3 (attenuation factor) of the downgoing 
ordinary and extraordinary waves. For details see text of figure 8. 




REAL PART OF EIGENVALUE 



I i I I . -- J. I 

0 0.2 0.4 0.6 0.8 1.0 

ANGLE OF INCIDENCE sin0 o 


Figure 10— Real part a of the eigenvalue of the upgoing acoustic-gravity waves versus 
the sine of the angle of incidence 6 q taking into account heat conduction, molecular 
viscosity, Coriolis force and ion drag. The parameter is the angular frequency 
(in sec - ^). The propagation conditions are: West to east propagation (full lines) and 
south to north propagation (dashed lines) at the equator. The dash-dotted lines sepa- 
rate different scales. 
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